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ABSTRACT 

The integrated brightness of the Sun shows variability on time-scales from minutes 
to decades. This variability is mainly caused by pressure mode oscillations, by granu- 
lation and by dark spots and bright faculae on the surface of the Sun. By analyzing 
the frequency spectrum of the integrated brightness we can obtain greater knowledge 
about these phenomena. It is shown how the frequency spectrum of the integrated 
brightness of the Sun in the frequency range from 100 to 3200 /iHz shows clear signs 
of both granulation, faculae and p-mode oscillations and that the measured character- 
istic time-scales and amplitudes of the acoustic signals from granulation and faculae 
are consistent with high-resolution observations of the solar surface. Using 13 years of 
observations of the Sun's integrated brightness from the VIRGO instrument on the 
SOHO satellite it is shown that the significance of the facular component varies with 
time and that it has a significance above 0.99 around half the time. Furthermore, an 
analysis of the temporal variability in the measured amplitudes of both the granu- 
lation, faculae and p-mode oscillation components in the frequency spectrum reveals 
that the amplitude of the p-mode oscillation component shows variability that follows 
the solar cycles, while the amplitudes of the granulation and facular components show 
signs of quasi-annual and quasi-biennial variability, respectively. 

Key words: Sun: faculae, plages - Sun: granulation - Sun: helioseismology - Sun: 
oscillations 



1 INTRODUCTION 

The frequency spectrum of observations of the integrated 
brightness of the Sun is often referre d to as the solar acous- 
' tic spectrum l|lben k, Ma haffv 1976), as it is mainly used to 
study pressure or sound waves inside the Sun. These sound 
waves can be observed on the surface of the Sun as oscil- 
lation pattens which can be represented as spherical har- 
monics of the eigenfunctions of the sound waves inside the 
Sun. These sound waves lead to a series of distinct peaks in 
the solar acoustic spectrum (see Fig. 1) and by studying the 
frequencies, amplitudes and line-widths of these peaks we 
can through helioseismology obtain kno wledge of the physic s 
that takes place inside the Sun (see e.g. iGough et al.lll996h . 
but more information are hidden in the solar acoustic spec- 
trum. At lower frequency the solar acoustic spectrum con- 
tains signals from spots, granulation and faculae. In fact the 
acoustic signal from granulation and faculae continues up 
to the frequencies of the p-mode oscillations and even to 
higher frequencies and therefore this signal is often referred 
to as the acoustic background. Studies of Sun-like stars have 
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shown t hat other stars show a similar acoustic background 
(see e.g. iMichel et af]|2008l : IChaplin et aill2010l ). 

lHarvevI (| 19851 ) could be the first to suggest that the 
solar acoustic background was modelled with a number of 
exponentially decaying functions. The idea is that the signal 
from e.g. granulation is understood as a random signal with 
some memory and thus its autocorrelation function should 
be given by an exponentially decaying functi on defined by 
a cha racteristic time-scale and an amplitude |Baudin et alj 
[2007]). Long uninterrupted observations from the South Pole 
and from the SOHO satellite have later revealed that the 
solar acoustic background contai ns signals from many more 
phenomena than just granulation. lHarvev et afl (|1993l) iden- 
tified a total of 6 difference components in the acoustic spec- 
trum from observations of the Ca n K line in the solar at- 
mosphere from the South Pole. The identified components 
were, from lower to higher frequency: non-periodic fluctua- 
tions attributed to active region evolution; non-periodic fluc- 
tuations attributed to granulation overshoot; non-periodic 
fluctuations attributed to chromospheric bright points; peri- 
odic fluctuations attributed to the photosphere (i.e. p-mode 
oscillations) and periodic fluctuations attributed to the chro- 
mosphere (i.e. chromospheric oscillations). Whereas the sig- 
nals from granulation and p-mode oscillations are reported 
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Figure 1. The solar acoustic spectrum. The spectrum was cal- 
culated as the mean of 119 one months spectra (see text). Both 
the granulation component (around 800 /^Hz), the facular com- 
ponent (around 1800 /^Hz) and the p-mode oscillations compo- 
nent (around 3000 /^Hz) are seen as changes in the slope of the 
spectrum. The different components are more clearly marked in 
Fig. 2. It is also seen that the spectrum declines faster at fre- 
quencies higher than the p-mode oscillations or the atmospheric 
acoustic cut-off frequency. The spectrum was calculated from 13 
years of observations of the Sun's integrated brightness from the 
VIRGO instrument on the SOHO satellite and the peak at 5555 
/iHz is an artefact. 



in all studies, the signals associate d with the chr o mosp here 
are only reported in the study by lHarvev et al. (1993), as 
this is the only study that analyzes the acoustic background 
in observations of the Ca II lines. The feature between the 
granulation signal and the p-mode oscillations has on the 
other hand also been reported in most studies, especially 
studies based on intensity measurements, but an agreement 
on t he origin of th i s feat ures has not been reached. 

lAigrain et all l|2004 ) analyzed observations of solar ir- 
radiance from the VIRGO instrument on SOHO and iden- 
tified not only a granulation component in the frequency 
spectrum, but also components from active regions, super- 
and meso-granulations and bright poi nts. The char a cteris - 
tic time-scale of the component that lAigrain et all (|2004l ) 
identified as bright points is si milar to the charact eristic 
time-scale of the component that lHarvev et all (Il993l) iden- 
tified as chromospheric bright points. Vazquez Ramio et al.l 
identified the same component in integrated light ob- 
servations from VIRGO, but speculated that this component 
might be due to a second granulation population - as shall 
be showen later this is not likely. 

Bright points and faculae are related phenomena and 
they are both related to the magnetic network on the so- 
lar surface. The difference is that whereas bright points 
are mainly seen in the inter-granular lanes, faculae are 
seen inside the granules and that whereas bright points 
can maintain their brightness over several tens of minutes 
l|Berger et al J 120071) . faculae chang e significantly on granu- 
lar evolution time-scales (|Pe Pontieu et aTl l2006t ). For this 
reason it is more likely that the component between the 
granulation and p-mode oscillation components is caused 
by faculae rather th an bright points. 

In the picture of lBerger et al.l l|200Tl ) faculae are deffined 
as "the edges of granules seen through the "forest" of the 



magnetic field in plages and the network" . The faculae does 
in other words occur in a localized region of the gran- 
ules and it is thus clear that the characteristic time-scale 
of faculae will be a bit shorter than that of granulation. 
There are no quantitative studies that directly relates the 
time-scale of granulation to the time-scale of faculae based 
on high-resolution observations. A qualitative argument for 
that the characteristic time-scale of faculae is shorter than 
th at of granulation is giv en in the small movie in Fig. 9 
of |Pe Pontieu et al.1 l|200fj) . This small movie of a 3D radia- 
tive magneto convective simulation shows 3 distinct granules 
that live throughout the full length of the movie (330 sec) 
and a number of faculae on these granules. If one look e.g. 
at the facular at (0'.'3,l'.'4) it is clearly visible at t = sec, 
but completely gone at t = 90 sec. At t — 120 sec a new fac- 
ular starts to evolve just to the right of where the other was 
located and at t — 210 sec, this facular has its peak bright- 
ness. In other words, though the evolution of faculae follows 
the evolution of the granules it is clear from the movie that 
the same granule can, during its lifetime, contain a num- 
ber of faculae - i.e. the characteristic time-scale of faculae 
is shorter than that of granulation. 

A few studies (|Trampedach et all Il998l ; iLudwigj 120061 ; 
iLudwig et al.l 120091 ) have used hydrodynamical models of 
stellar atmospheres to simulate acoustic spectra of Sun-like 
stars, but so far no such studies have been made using mag- 
netohydrodynamical models, which is needed in order to re- 
produce the facular component. 

In this paper we will use observations from the 
green channel on the VIRGO instrument on SOHO 
( Froh lich et al.l Il997l ) to analyze the acoustic signature of 
faculae. We will show how the characteristic time-scale and 
amplitude that we measure of this component in the acous- 
tic background are consistent with measurements based on 
high-resolution observations of the solar surface. We also an- 
alyze any temporal variations in the different components in 
the solar acoustic background - including the facular com- 
ponent. 

The paper is arranged as follows. In Section 2 we give 
a short theoretical description of the acoustic background 
and present the model that we use to model the granula- 
tion, the faculae and the p-mode oscillations components in 
the acoustic spectrum. This model is used in the analysis 
in Section 3, where we also describe a statistical test of the 
significant of the facular component in the acoustic back- 
ground. Section 4 presents the results of the analysis and 
concluding remarks are found in section 5. The paper also 
includes an appendix that describe a test based on Monte 
Carlo simulations of the statistical method developed in Sec- 
tion 3. 



2 THEORY 

T he original form ulation of the solar acoustic background 
by lHarvevI (Il985l ) is the following: 



(i) 



1 + {2-kvt) 2 ' 

where f(v) is the power density at frequency v, a is the am- 
plitude of the signal and r is the characteristic time-scale 
Eq. 1 is normalized to Parseval's theorem so that the vari- 
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Figure 2. The same a Fig. 1, but here with the model containing 
a facular component overlaid. The solid line shows the full model 
including all three components, whereas the dashed lines show 
the different components in the model (i.e. granulation, faculae 
and p-mode oscillations). The parameter values used in the model 
are given in Table 1. 



ance of the signal in the time domain equals the sum of f(v) 
over all frequencies: 



(2) 



which means that a in eq. 1 does in f act reflects the var i- 
ance of the signal in the time domain (|Baudin et al.|[2007t ). 
Note that every component in the acoustic background will 
contribute with the signal given by eq. 1 to the spectrum. 
The spectrum will thus be a sum of individual signals from 
i.e. granulation and faculae each defined by unique values of 

a and r. 

Unfortunately, the original model bv lHarvevI (1 19851 ) fails 
to reproduce the observed solar acoustic background for fre- 
quencies higher than the atmospheric acoustic cut-off fre- 
quency. The reason for this is th at granulation cannot be 
modelled with turbulent cascades (Nordlun d et al.lll997l ) as 
it is done in the drift model by Irlarvevl (| 19851 ). Turbulence 
shows a distribution with a slope of around —2 in power, 
convection on the other hand has a lower limit in the time 
domain on which changes can take place. This means that 
on small time-scales (or at high frequency) convection is not 
noisy whereas turbulence is and therefore the acoustic back- 
ground decays with a slope smaller than —2. In the Sun the 
decay rate turns out to be around —4 at frequencies higher 
than the atmospheric acoustic cut-off frequency. Taking this 
into account the Harvey model extends to: 



/(") = 



(3) 



1 + (2-kvt) 2 + {2-kvtY' 

iKarofj I 20081 ) used this model to successful describe the 
granulation and facular components in the solar acoustic 
background between 400 ^iHz and up to 8000 /iHz, which is 
clos e to the Nyqu i st freq uency. 

lHarvev et alJ l|l993l ) used another approach to solve the 
problem with the fast decline at high frequency and consid- 
ered the slope as a free parameter. When this is done the 
decline will often be faster then —2 a nd this will so l ve the 
problems at high frequency. In this wav lHarvev et all ij 19931 ) 
found slopes of —5.6 and —5.0 for the granulation and the 



Figure 3. The same as Fig. 2, but here the model does not 
contain a facular component. The parameter values used in the 
model are given in Table 1. 



bright point components, respectively. The slope gives the 
decay rate, which 'calibrates' the amount of memory i n the 
process responsible for the noise (i Harvev et al.1 1993) and 
there is no reason that this value should be exactly —2, or 
that is should be the same for the different components. On 
the other hand, there is a physical argument why the slope 
should not be the same at low and high frequency (which 
were given above), so allowing the slope to be a free param- 
eter only partly solves the problem. 

We have therefore chosen another approach and mod- 
elled only the lower frequency part of the spectrum (up to 
3200 /iHz). This means that no extra component is needed 
in order to model the high-frequency part. The following 
model is thus used to model the granulation and facular sig- 
nals in the solar acoustic background between 100 and 3200 
/iHz: 



/O) 



1 + (2-kvt) 



(4) 



The envelope of the p-mode oscillations shows variabil- 
ity on a weekly time-scale with an amplitude of 6.2%, which 
is properly mainly due to changes in the d amping rates of 
the oscillation modes (Kield sen et all 120081 ) . The envelope 
of the p-mode oscillations therefore needs to be modelled 
along with the acoustic background in order to disentangle 
changes in the damping rates of the oscillation modes from 
changes in the acoustic background. A Gaussian envelope is 
therefore also included in the model in order to model the en- 
velope of the p-mode oscillations, but as the high frequency 
part of the spectrum is not analyzed it is not necessary to 
include a white noise background. Of course the background 
will have some effect on the amplitudes that are measured 
of the granulation and facular components, but as the white 
noise level in the VIRGO observations is generally around 
three orders of magnitudes smaller than the granulation and 
facular amplitudes this effect is negligible. 

A Gaussian has here been chosen for the envelope of the 
p-mode oscillation s, but one or two Lorentzians could als o 
have been chosen (|Harvev et al ] |l993l : iLefebvre et aT]|2008l ) . 
This would how ever not affect the final results significantly 
and it is not clear which functional form that agrees best 
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with the behavior of the Sun and other Sun-like star s in 
general (see e.g. the simulations bv lHoudek et al.1l 19991 ). 

Figs. 2 & 3 show examples of fits to the solar acoustic 
background using the models described above either con- 
taining or not containing a facular component. The obser- 
vations and the calculation of the spectrum are described in 
Section 4. 



3 ANALYSIS 

The data analysis consists of three steps; calculating the 
power density spectra, modelling the power density spec- 
tra and evaluating the significance of the different measured 
parameters and the goodness of the different models. 



3.1 Observations 

For the analysis we use a 13-year time series from the 
green channel on the VIRGO instrument of SOHO l|Frohlicbl 
l2009h . There are some gaps in this time series, especially 
around the SOHO 'vacation' and recent measurements do 
suffer from increased noise level. In the analysis performed 
here we have divided the full 13-year time series into sub- 
strings of one month each. In order to get the lowest prob- 
ability for artefacts such as gaps and increased noise levels 
to affect the data analysis we have performed a discrimina- 
tion of both bad data points and bad substrings. Bad data 
points were identified by calculating a logarithmic running 
box variance log (<7i): 



log (cTi) 



loe 



l 

N 



j — i-\-w 

E i 



(5) 



where Xj is the jth observation, fi is the mean value of the 
observations between i — w and i + w and w is the width of 
the running box which was set equal to 10 points (N = 21). 
Using this formulation all data points with a logarithmic 
running box variance that deviated more than 5a from the 
mean were then removed. This removed around 5% of the 
data points. The reason for using the logarithmic variance 
is that many data points show artificially low scatter. 

After removing the bad data points all monthly sub- 
strings with a duty cycle lower than 95% were also removed. 
This left us with 119 monthly substrings out of 147 possi- 
ble. We did test that the relative high demands for selecting 
good data points did not affect the results. The test showed 
slightly divagating results when using either the 119 good 
monthly substrings or all the 147 substrings, but the differ- 
ences were insignificant. 

The power dens i ty spectra w ere calculated using least- 
squares (|Lomblll976l ; K aroffeoolj ). Each spectrum was nor- 
malized by the effective observation length given as the re- 
ciprocal of the area under the window function in order to 
convert the spectra into power density spectra. 



3.2 Modelling the acoustic background 

Acoustic backgrounds have earlier been modelled by 
smoothing the spectra with a Gaussian running mean, with 
a width of a few hundred /iHz, and then assuming that the 



differences between the smoothed observed spectrum and 
the model followed a normal distribution, which allowed a 
comparison utilising least squares (se e e.g. iKieldsen et alj 
l200a lKarof3l2008l ; IChaplin et alj|20ld) . Though this works 
in practice, there are no mathematically formal arguments 
for that the differences between the smoothed observed spec- 
trum and the model should follow a normal distribution. A 
method for modelling the observed spectra that is statis- 
tically valid and that does not require the spectra to be 
smoothed with a Gaussian running mean is therefore pre- 
sented below. 

In order to model the solar acoustic background we 
make use of maxim um likelihood estimators (see e.g. 
I Anderson et al.lll990l . and references herein). The purpose 
is to calculate a logarithmic likelihood function £ between 
iV independent measurements Xi - i.e. power density at a 
given frequency, and a model fi given by a set of param- 
eters A. A dopting the notation of iToutain fc Appourchauxl 
(| 19941 ) and lAppourchaux et al.l (1 19981 ) the logarithmic like- 
lihood function can be calculated as: 



(6) 



where p(xi, A) is the probability density function. The posi- 
tion of the minim um of £ in A space will g ive us the most 
likely value of A l|Appourchaux et al.l K1998) . The so-called 
formal error bars can be calculated as the diagonal elements 
of the inverse of the Hessian matrix h: 



d 2 £ 

dXidXj 



(7) 



In order to calcula te the probabil i ty den sity function we 
use the formulation by I Appourchauxl (|2004l ) where a binned 
version of the origi nal spectrum is com pared to the model. 
The formulation bv l Appourchauxl l|2004l ) can be used on any 
model fi and not just for calculating the significance of peaks 
in the spectrum. We thus calculate the binned spectrum 
Si(x,n) by summing Xi over n bins: 



k = i+n/2 



5l (x,n) = -l T . Yl 



(8) 



k=i—n/2 



Following lAppourcha ux (2004) the probability density func- 
tion is then given as: 



p[Si(x,n),fi(X)} = —^-S t (x,n) 
1 \yij 



Vi-\ Si(x,n) 



(9) 



where T is the Gamma function and [ii and Vi are given as: 



Hi = 



and 



fc = i + n/2 

(n + 1) f"W 

k=i — n/2 
k = i+n/2 ' 

£ /2(a) 

k=i — n/2 



(10) 
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Table 1. Acoustic background parameters 





faculae 


no faculae 


CT gran 


73.0 ± 1.6 ppm 


85.8 ± 1.5 ppm 


Tgran 


214.3 ± 2.9 sec 


219.0 ± 4.9 sec 


a 


-3.5 ± 0.3 


-1.8 ± 0.1 


°"fac 


61.5 ± 1.9 ppm 




Tfac 


65.8 ± 2.5 sec 




a 


-6.2 ± 0.7 




Hose 


0.7 ± 0.04 ppm 2 //xHz 


0.52 ± 0.05ppm 2 //xHz 




3104 ± 36 fiHz 


3084 ± 44 fiHz 


Width 


316 ± 36 fiHz 


254 ± 32 /iHz 



k = i+n/2 
k=i — n / 2 



k=i+n/2 
k=i — n/2 



(11) 



n appears in p and v in order to get the right amplitudes. 

3.3 Testing the facular hypothesis 

Using the formalism above the model in eq. 4, either contain- 
ing or not containing a facular component, can be compared 
to an observed spectrum. The model containing a facular 
component is based on a total of 9 free parameters - i.e. 
three from granulation, three from faculae and three from 
the p-mode oscillations, whereas the model not containing a 
facular component is based on only 6 free parameters. Hav- 
ing the maximum likelihood of each model from eq. 6 we 
are then able to calculate the statistical significance of the 
facular component in the acoustic spectrum by calculating 
the logarithmic likelihood ratio A: 

In A = £(\ p+q )-£(\ p ). (12) 

In our c ase p equals 6 and q equals 3. Following I Wilksl (I1938T ) 
(see also lAppourchaux et al.lll99cf ) we can then compare the 
value of — 21nA to a x 2 distribution with q degrees of free- 
dom and in this way calculate the significance of the facular 
component in the acoustic spectrum. 



4 RESULTS 

Two kinds of results are presented: One set of results on the 
average spectrum of the 119 one months substrings and one 
set of results on the temporal variability in the spectra of the 
individual substrings. The analysis is made by binning 100 
frequency bins (n = 100), which gives a frequency resolution 
in the binned spectra of around 39 //Hz. Error bars for all 
parameters were calculated as the diagonal elements of the 
inverse of the Hessian matrix (see eq. 7) 



4.1 Average results 

The average spectrum from the 119 one month substrings is 
modelled with a model containing a facular component and a 
model not containing a facular component. The results from 
the modelling are given in Table 1 and the average observed 



spectra with the models overlaid are shown in Figs. 2 & 3. 
As the decay rates were treated as free parameters in both 
models the decay rates for the granulation component are 
different in the two models. This leads to differences in the 
models at high frequency as the obtained decay rate a of 1.8 
±0.1 for the granulation component in the model without 
a facular component is too small to describe the acoustic 
background at high frequency. 



4.2 Comparison to high-resolution observations 

A number of studies have used high-resolution observations 
of the solar surface to measure the amplitude and time-scale 
of granulation and faculae, but in general there is little agree- 
ment between the different results. The reason for this is 
probably differences in the way the amplitude and the time- 
scale are measured, but also differences in the resolution of 
the observations, differences in the solar activity level at the 
time of the observations and whether the observations were 
taken close to the disk centre or close to the limb. The high- 
resolution observations can therefore only be used as a rough 
comparison to the seismic results. 

If a typical size o f granules (or cell area) of 1.2 
Mm 2 l|Del Morol l2004h is assumed, then around 2.5 
million granules are found on the visible solar sur- 
face. The contrast of a typical granule is aroun d 20% 
l|Wedemever-Bohm fc Rouppe van der Voortl I2009T ) . The 
typical contrast of faculae is around 50 % at the peak dis - 
tance from the disk center (fj, ~ 0.6) dBerger et alj 120071 ) 
and the relat ive fraction of gra nules that show faculae is 
around 1.8% l|Berger et al.lll995f ). 

If each granule on the solar surface is assumed to cause 
a signal of relative amplitude ^ at a random time and the 
entire solar surface is covered with granules, then this will 
ca use a time ser i es of integrated Sun light with a variance 
of (iLudwid 1200(1 iLudwig et ai]|2009l ,: 



granule 



yj -^granule 



(13) 



where N granu i e is the number of granules present on the 
visible solar surface at any given time. 

The amplitudes of the granulation component were 
measured to be 73.0 ± 1.6 and 85.8 ± 1.5 ppm for the model 
containing a facular component and the model not contain- 
ing a facular component, respectively. This is consistent with 
the theoretical estimate of around 126 ppm that is found by 
dividing 20% with the square root of 2.5 million granules. 

The amplitude of the intensity signal from faculae fol- 
lows that of granules, except that it cannot be assumed that 
the entire solar surface is covered with faculae. In the most 
simple approach the amplitude should therefore be scaled 
with the relative coverage of faculae a " : 



( SL ) 

Nfa 



granule N granu l e 



(14) 



The amplitude of the facular component was measured 
to be 61.5 ± 1.9 ppm. This agrees roughly with the theo- 
retic estimate of around 42 ppm obtained by multiplying 
~50% with the square root of ~1.8% and dividing it with 
the square root of 2.5 million granules. 
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Figure 4. The logarithmic likelihood ration between the model 
with and without the facular component as a function of time. No 
significant variability following the solar cycle is seen. The solid 
black line shows a Gaussian running mean with a width (a not 
FWHM) of 0.4 year. 



Figure 5. Temporal variability of the sunspot number for com- 
parison. The figure shows the international sunspot number ob- 
tained from the Solar Influences Data Analysis Center at the 
Royal Observatory of Belgium. The solid black line shows a Gaus- 
sian running mean with a width of 0.4 year, as in the other figures. 



The characteristic time-scales of convection were mea- 
sured to be 214.3 ± 2.9 and 219.0 ± 4.9 sec for the model 
containing a facular component and the model not contain- 
ing a facular component, respectively. This is consistent with 
the coherence time-scale between 156 and 246 sec found by 
Del Moro (2004). For faculae a characteristic time-scale of 
65.8 ±2.5 sec was measured which is consistent with the fact 
that it should be shorter than the granulation time-scale. 

The amplitude and characteristic time-scale of the gran- 
ulation and facular components in the acoustic spectra are 
consistent with the high-resolution measurements. This in- 
dicates that the interpretation that the second bump in the 
acoustic spectrum originates from faculae is consistent - 
though the comparison does not make the argument alone. 
The interpretation is also based on the relative relations be- 
tween the amplitude and characteristic time-scale of granu- 
lation, faculae and p-mode oscillations. 

Eq. 13 shows why the component in the acoustic spec- 
trum between the granulation and p-mode oscillation com- 
ponents does not origin ate from a second granulatio n pop- 
ulation as suggested by I Vazquez Rarnio et adj J2005 ). This 
suggestion is based on the work by iDel Morol |2004h who 
finds, using high-resolution images of quiet granulation, a 
second population of granules with a shorter coherence time, 
but also smaller brightness contrasts. The fact that this sec- 
ond population of granules is both more numerous and have 
smaller brightness contrasts than ordinary granules means 
that it should not be visible in the acoustic spectrum as the 
numerator of eq. 13 would be smaller and the denominator 
would be larger than for ordinary granulation. 



4.3 Statistical significance of the facular signal 

The significance of the facular component in the average ob- 
served spectrum can be calculated from the likelihood ratio 
between the two different models. The significance was mea- 
sured to be 83%. This means that though the model with the 
facular component and a total of 9 free parameters does de- 
scribe the average observed spectrum better than the model 
without the facular component and 6 free parameters, it 



does not describe it significantly better. Visually, it is clear 
that the fit in Fig. 2 is better than the fit in Fig. 3, but given 
the statistic of the observed spectra and the number of free 
parameters it is not significant better. 

In order to analyze why the facular component is not 
significant in the average observed spectrum we calculate 
the significance of the facular component in the spectrum 
of each individual one month substring. The decay rates 
and characteristic time-scales for the granulation and fac- 
ular components and the frequency of maximum power for 
the p-mode oscillation component had to be fixed to the val- 
ues given in the first column of Table 1. If this was not done 
the location and slope of the exponential decaying function 
would change in order to make up for changes in ampli- 
tude. This means that the models used for modelling the 
119 individual one month substrings contained 3 and 4 free 
parameters, respectively. 

The logarithmic likelihood ratios between the models 
with and without a facular component for the 119 individ- 
ual one month substrings are shown in Fig. 4. It is seen 
that most of the logarithmic likelihood ratios are between 
and 10 with a few scatter points below and above. The 
figure does reveal a few points with logarithmic likelihood 
ratios below 0. A negative logarithmic likelihood ratio means 
that in these months the model without the facular compo- 
nent describes the observed spectrum better than the model 
with the facular component, or in other words, that adding 
more free parameters to the model makes a worse agree- 
ment between the observed spectrum and the model for 
these months. The reason for this is that the decay rates 
of the granulation components are not the same in the two 
different models and for some months the slow decay rate in 
the model without the facular component describes the ob- 
served spectrum better than the fast decay rate in the model 
with the facular component. These scatter points with log- 
arithmic likelihood ratios below could be the reason for 
the facular component not being significant in the average 
observed spectrum. 

Clearly, a lot of variability is seen in Fig. 4, but no 
trend known to follow the visibility or amplitude of faculae 
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Figure 6. Histogram of the measured significances of the facular 
component in the 119 one months substrings. The significances 
are calculated using the logarithmic likelihood ratios in Fig. 4 
and assuming that these ratios follow a x 2 distribution with one 
degree of freedom. It is seen that 56 months (47%) have a signifi- 
cances higher than 0.99, while 75 (63%) have a significance higher 
than 0.95. 
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is seen. In particular no clear variability that follows the 
solar cycle is seen (see Fig. 5 for comparison). The reason 
for this could either be that the measured ratios are not 
accurate enough to show intrinsic facular variability or that 
the acoustic signal from faculae depends stronger on the 
number of available granules than on the fraction or contrast 
of faculae (see eq. 14). 

Fig. 6 shows the distribution of the significances of the 
facular signal measured in the 119 one month substrings 
assuming that the logarithmic likelihood ratios follow a \ 2 
distribution with one degree of freedom. This distribution 
shows a clear increase towards a significance of one that falls 
of rapidly down to a significance of around 0.97 and then 
have a long tail down to smaller significances (also smaller 
significances than what is shown on the figure) . The conclu- 
sion from the figure is that the acoustic signal from faculae 
is significantly present in the observed solar acoustic back- 
ground 47% of the time. 

This figure support our impression that the acoustic sig- 
nal from faculae is present on the Sun most of the time and 
that the reason that the average spectrum does not shows a 
significant facular signal is because of a few months where 
the significance is really low - i.e. where the model with- 
out the facular component describes the observed spectrum 
better. The same effect is seen in the spectrum calculated 
from the full 13 years of observations - i.e. that the acoustic 
signal from faculae in not significantly present in the spec- 
trum. Again, the reason for this is most likely that the signal 
is not present all the time, which reduces the significance of 
the signal when analyzing the full 13 years of observations. 



4.4 Temporal variation of p-mode oscillations 

Both the frequencies, amplitudes and line-widths of the 
individual p-mode oscillation modes in the acoustic spec- 
trum are known to change with the solar cycle (see e.g.: 
iLibbrecht fc Woodardl Il990l : IChaplin et all l200Gh and we 
therefore also expect the average shape and position of the 



Figure 7. Temporal variability of the height (left) and width 
(right) of the p-mode oscillation envelope for the model with a 
facular component (top) and the model without a facular com- 
ponent (bottom). The solid black lines show a Gaussian running 
mean with a width of 0.4 year. As expected a clear anti-correlation 
is seen between the measured height of the p-mode oscillation en- 
velope and the solar cycle. 



p-mode oscillation envelope to change with the solar cycle. 
The temporal variability of the height and width of the p- 
mode oscillations envelope is shown in Fig. 7 for the model 
with and without a facular component (the position of the 
envelope is kept fixed in both models). There is a good agree- 
ment between the results from the two different models - 
with a linear Pearson correlation coefficient of 0.95 (0.97 
for the smoothed curves) for the height and 0.83 (0.70 for 
the smoothed curve) for the width. A clear anti-correlation 
is seen between the height of the envelope and the solar 
cycle (see Fig. 5 for comparison) so that the height has min- 
imum around 2000, which coincides with the peak of the 
solar cycle. The anti-correlation between the height of the 
envelop e and the solar cycles is consistent with other stud- 
ies (i.e. iKieldsen et a l. 2008; Garcia et al.ll201ol ) and can be 
explained by the fact that the life-time of the p-mode os- 
cillation is known to decrease with increasing solar activity 
jChaplin et alj l2003h . The envelope width does not show 
any correlation with the solar cycle. The position of the en- 
velope, which is not analyzed here, as it was fixed to the 
average value, does only show small ch anges compared to 
the hight and width (|Garci'a et alj|2010l ). 

The uncertainties on the parameters measured in the 
average spectrum are comparable in size to the uncertain- 
ties measured in the individual one month substrings. This 
reflects that the scatter seen in i.e. Figs. 7 & 8 is much larger 
than the uncertainties on the individual measurements - in 
other words it reflects that the parameters does show intrin- 
sic variability on a monthly time-scale. 
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Figure 8. Temporal variability of the amplitude of the granu- 
lation component for the model with a facular component (top) 
and without a facular component (bottom). Again, the solid black 
lines show a Gaussian running mean with a width of 0.4 year. 
Some signs of a quasi-annual variability is seen, especially in the 
results from the model with the facular comp onent, as was also 
the case in the study by iLefebvre et alj {2008 



4.5 Temporal variation of granulation 

The temporal variability of the amplitude of the granula- 
tion component does not show any clear correlation with 
the solar cycle, neither for the model with or without the 
facular component. The amplitudes of the granulation com- 
ponent show an offset between the two different models and 
only some agreement is seen in the temporal variability - 
with a linear Pearson correlation coefficient of 0.74 (0.84 for 
the smoothed curves). Instead some signs of a quasi-annual 
modulation is seen, especially in the results from the model 
with the facular component. The ratio between the variance 
of the measurements (the scatter) and the mean value of 
the uncertainties is 5.12 for the model with and 4.97 for the 
model without the faculae component, indicating that the 
variability is not due to noise. A quasi-annual modulation 
of both the granulation characteristic time-scale and ampli- 
tude was also seen in the analysis of observations from the 
GOLF instrument on SOHO by ILefebvre etafl ll2008h in 
phase with the variability seen here. ILefebvre et all (|2008t ) 
explained the modulation by the fact that the observation 
altitude (or depth) of the GOLF instrument in the solar 
photosphere changes with the observation wavelength which 
again changes with the orbital velocity of SOHO, which 
changes with a period of one year. As the VIRGO obser- 
vations analyzed here are intensity observations, the effect 
of a wavelength change should not be prominent, but the 
quasi-annual variability is still seen, especially around solar 
activity maximum. 



uncertainties are relatively large for these measurements, 
and the measured amplitudes show only little scatter within 
these uncertainties. The amplitude of the scatter is around a 
few ppm. The Gaussian running mean of the measurements 
(represented with the black line) show a quasi-biennial vari- 
ability with maximum around 2001, 2003, 2005 and 2007. 
The amplitude of this variability is around one ppm - i.e. 
significantly lower than the scatter of the individual points. 
It is therefore not clear if the quasi-biennial variability is 
intrinsic to the Sun or if it is just random noise. 

A similar quasi-biennial variability was seen in the resid- 
ua ls of the p-mode oscil lations frequency shifts in the study 
bv lFletcher et al.ll|2010D . In Fig. 10 we therefore compare the 
residuals of the measured amplitudes of the facular compo- 
nent with the residuals of the p -mode oscillations frequency 
shifts from lFletcher et"aL l (|2010l ). In order to do that we have 
subtracted a 3rd order polynomial from the measurements 
in Fig. 9 and binned the one month measurements to the 
temporal resolution of the residuals of the p-mode oscilla- 
tions frequency shifts (182.5 days). The two sets of measure- 
ments show a correlation with a linear Pearson correlation 
coefficient of 0.34. The reason that no strong correlation is 
seen could however be the relative large uncertainties on 
both sets of measurements and it can thus not been ruled 
out that the faculae and the p-mode oscillations frequency 
shif ts have the same orig in . 

iFletcher et al.l (|2010l ) speculate that the quasi-biennial 
variability seen in the residuals of the p-mode oscillations 
frequency shifts originates from a second solar dynamo lo- 
cated not at the so-called tachocline at the base of the 
convection zone, but in the strong shear layer just below 
the surface. This idea of a near surfa c e dyn amo, which has 
also been suggested by iBrandenburel ((2005 ) , could also ex- 
plain the short periodic magnetic var iability that has re- 
cently observed in other Sun-like stars |Bohm-Vitensell2007l ; 
iGarcia et all |2010| ; iMetcalfe et~aH |2010| ). If the presents of 
faculae and p-mode oscillation frequency shifts are indeed 
correlated (which we can however not claimed based on the 
present measurements) then this could suggests that the 
presence of faculae on the surface of the Sun is controlled 
not only by the 11-years solar cycles, which is expected to be 
located in the solar tachocline, but also by a quasi-biennial 
dynamo located i n the strong radial shear layer just below 
the solar surface (|Brandenburdl2005l ) 

In general it has to be noted that the quasi-annual and 
quasi-biennial periods, seen in the granulation and facular 
amplitudes, are very peculiar. Varying instrumental effects 
such as altitude, pointing, temperature or the like could be 
the culprit. One of the things that were done in order to 
test this was to redo the analysis using only the monthly 
substrings where the facular signal had a significance larger 
than 95 and 99% and this did not change any of the conclu- 
sions on the facular signal - i.e. the same kind of periodicity 
were still seen. 



4.6 Temporal variation of faculae 

Though the facular component is not significantly present 
in the acoustic spectra of all the one month substrings the 
amplitude of the component can still be measured in all 
the substrings. These amplitudes are plotted in Fig. 9. The 



5 CONCLUSION 

13 years of observations of the integrated brightness of the 
Sun from the VIRGO instrument on SOHO have been an- 
alyzed. Three distinct features are seen in the frequency 
spectrum of these observations. By comparing the measured 
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Figure 9. The temporal variability of the amplitude of the fac- 
ular component. Again, the solid black lines show a Gaussian 
running mean with a width of 0.4 year. Though the measured am- 
plitudes show only little scatter within the uncertainties the black 
line reveals a quasi-biennial variability. A similar quasi-biennial 
variability is also seen in the res iduals of the p-mode oscillations 
frequency shifts in the study by [Fletcher ct al. (2010) 

, which suggest that the variability in the amplitude of the facular 
component could be caused by a second dynamo. 

characteristic time-scale and amplitude of these phenom- 
ena to high-resolution observations of the solar surface it is 
shown that the three features can convincingly be explained 
as granulation, faculae and p-mode oscillations. 

The feature in the acoustic background related to fac- 
ulae is the weakest and it is not significantly present in the 
average spectrum made from 119 individual one month spec- 
tra. The feature is on the other hand present in 56 (47%) 
out of the 119 one month spectra with a significance larger 
than 0.99 and in 75 (63%) with a significance larger than 
0.95. 

The temporal variability of the amplitudes of the gran- 
ulation, facular and p-mode oscillations components in the 
solar acoustic background have also been analyzed. The p- 
mode oscillations show a clear anti-correlation with the solar 
cycle and the amplitude of the granulation component show 
some signs of an quasi-a nnual variability, whic h is also seen 
in velocity observations (Lef ebvre et alj [20081 ). but no cor- 
relation with the solar cycle. The amplitude of the facular 
component does also not show any correlation with the solar 
cycle, but instead some signs of a quasi-biennial variability 
is seen. The residuals of the temporal variability of the am- 
plitude of the facular component show a correlation with a 
linear Pearson correlation coefficient of 0.34 with the residu- 
als ofJJiej3-m£de osdllations frequency shifts from the study 
bv lFletcher et all (|2010l ). 
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APPENDIX A: MONTE CARLO SIMULATIONS 
OF THE ACOUSTIC SIGNAL FROM FACULAE 

In order to test the robustness of the statical analysis devel- 
oped in Section 3, a number of Monte Carlo simulations were 
made. In the simulations it was assumed that the acoustic 
spectrum followed the model given in eq. 4 and an acous- 
tic spectrum was simulated for the parameters given in Ta- 
ble 1. Using an inverse Fourier transform the acoustic spec- 
trum was convented from the frequency to the time domain, 
where realistic random noise were added at each time-step. 
Hereafter the acoustic spectrum was convented back into 
the frequency domain, where it was analyzed - the back- 
ground parameters and the significance of the facular signal 
measured. 

The two parameters that were changed in the simulation 
were the amplitude of the facular component and the noise 
level - all other parameters were fixed to the values given in 
Table 1. This means that for each set of facular amplitude 
and noise level the simulation returned not only measured 
facular amplitude and noise level, but also measured values 
for all the other parameters in Table 1, which could be used 
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Figure Al. Relation between the measured amplitude and the 
measured significance of the facular component in the simulations 
(here just given as the logarithmic likelihood ratio in order to en- 
hance the structures in the plot). It is seen that the relation is far 
from linear. The reason for this is, that the likelihood between 
the simulations and the model that does not include a facular 
component is not a linear function of the amplitude of the facu- 
lar component. The likelihood between the simulations and the 
model that does include a facular component is on the other hand 
linear. 

to evaluate the precision of the procedure. The amplitude 
of the facular component was changed between 45 and 70 
ppm and the noise level from to 80 ppm per measurement 
(equivalent to the point-to-point scatter in the time series) . 

The main result of this Monte Carlo simulations is 
shown in Fig. Al, which shows the relation between the am- 
plitude of the facular component and the logarithmic likeli- 
hood ration between the model with and without the facular 
component. It is seen that there is no linear correlation be- 
tween the two. This is a bit surprising as a spectrum with 
a hight amplitude facular component is expected to have 
a higher likelihood with a model with a facular component 
than a spectrum with a low amplitude facular component 
would have. This is also the case, but the problem is that 
there is no linear relation between the likelihood and the am- 
plitude of the facular component for the model that does not 
include the facular component. In other words, is the reason 
why no linear correlation is seen that though the likelihood 
between the model and the observations increases linear as a 
function of facular amplitude for the model with the facular 
component, the same (or more correctly the opposite) is not 
true for the model without the facular component. 

This does first of all explain why some sings of variabil- 
ity can be seen in Fig. 9, which shows the measured facular 
amplitude as a function of time, while not in Fig. 4, which 
shows the likelihood ration as a function of time. 

It does also explain why we can measure the amplitude 
and characteristic time-scale of the facular component and 
the uncertainties on these parameters even in months where 
the component is not significant. The reason is that in order 
to measure the amplitude and characteristic time-scale we 
only use the model that includes the facular component, 
while in order to measure the likelihood ratio we also uses 
the model that does not include the facular component. 

No trends were seen between the measured parameters 
and the noise level in the simulations, except for increased 
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Figure A2. Distribution of the difference between the input and 
output of the simulations. The figure shows a histogram of the 
difference between the facular amplitude that was used in the sim- 
ulations and the amplitudes that were returned from the model, 
in therms of the uncertainty returned from the modeling of the 
simulated spectra (solid line). The dashed line shows a Gaussian 
with a variance of one for comparison. The fact that the distri- 
bution is not to different from a Gaussian with a variance of one, 
reflects that the uncertainties returned from the modelling are 
realistic. 

scatter. For some of the simulations with the highest noise 
levels it was not possible to get the model to convert to the 
observed spectrum. 

The Monte Carlo simulations have also been used to test 
the uncertainties on the different parameters - especially the 
parameters related to the facular component. This is done by 
plotting a histogram of the difference between the parameter 
values used in the simulation and the measured parameter 
values divided by the uncertainty on the measured parame- 
ter values ( x °bs~ x tru C y jr t ne errors on the parameters were 
random distributed around a mean value we would expect 
that these histograms would all be Gaussian function with 
a variance of one. This is also the case for the parameters 
in Table 1 (the parameters related to the facular compo- 
nent, is shown in Fig. A2), though the variance might be 
a bit larger than one, this is not unexpected as the uncer- 
tainties are only formal uncertainties. In order to get more 
reliable uncertainties one would eventual have to preform a 
full Bayesian analysis as it is sometimes done when model- 
ing the p-modes using Marko v chain Monte Carlo (see e.g. 
lHandberg fc Campantell201ll 'l. 



